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Abstract. We present a neutron triple- axis and resonant spin-echo spectroscopy 
study of the spin correlations in untwinned YBa2Cu306+cc single crystals with x = 0.3, 
0.35, and 0.45 as a function of temperature and magnetic field. As the temperature 
T ^ 0, all samples exhibit static incommensurate magnetic order with propagation 
vector along the a-direction in the Cu02 planes. The incommensurability 6 increases 
monotonically with hole concentration, as it does in La2-xSra;Cu04. However, 5 
is generally smaller than in La2-2;Sr2;Cu04 at the same doping level, and there is 
no sign of a reorientation of the magnetic propagation vector at the lowest doping 
levels. The intensity of the incommensurate Bragg reflections increases linearly with 
magnetic field for YBa2Cu3 06.45 (superconducting = 35 K), whereas it is field- 
independent for YBa2Cu306.35 (Tc — 10 K). These results fit well into a picture in 
which superconducting and spin-density-wave order parameters coexist, and their ratio 
is controlled by the magnetic field. They also suggest that YBa2Cu306+x samples with 
X ^ 0.5 exhibit incommensurate magnetic order in the high fields used for the recent 
quantum oscillation experiments on this system, which likely induces a reconstruction 
of the Fermi surface. We present neutron resonant spin-echo measurements (with 
energy resolution ^ 1 /xeV) for T ^ that demonstrate a continuous thermal 
broadening of the incommensurate magnetic Bragg reflections into a quasielastic 
peak centered at excitation energy E = 0, consistent with the zero-temperature 
transition expected for a two-dimensional spin system with full spin-rotation symmetry. 
Measurements on YBa2Cu3 06.45 with a conventional triple-axis spectrometer (with 
energy resolution ~ 100 /ueV) yield a characteristic crossover temperature Tsdw ~ 30 
K for the onset of quasi-static magnetic order. Upon further heating, the wavevector 
characterizing low-energy spin excitations progressively approaches the commensurate 
antiferromagnetic wave vector, and the incommensurability vanishes in an order- 
parameter-like fashion at an "electronic liquid-crystal" onset temperature T^hc 
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150 K. Both TsDW and Telc increase continuously as the Mott-insulating phase is 
approached with decreasing doping level. These findings are discussed in the context 
of current models of the interplay between magnetism and superconductivity in the 
cuprates. 
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1. Introduction 

The phase diagram of the cuprate high-temperature superconductors includes three 
generic thermodynamic phases at zero temperature: a Mott-insulating phase with 
commensurate antiferromagnetic order, a d-wave superconducting phase, and a metaUic 
phase at low, intermediate, and high doping levels, respectively. In the limits of zero 
and large doping levels, the low-energy excitations of these phases are well understood. 
In the Mott-insulator, spin and charge excitations are separated by a large optical gap, 
and neutron scattering experiments have shown that the low-energy spin excitations are 
well described by a Heisenberg spin Hamiltonian [H In the metallic phase at high 
doping levels, various transport experiments indicate that the Fermi liquid theory yields 
an adequate description of the coupling between spin and charge excitations 0, IH E] • 
In the (i-wave superconducting (SC) phase, however, a quantitative description of the 
coupling between spin and charge excitations remains a challenge for current research, 
and there are some indications that such a description will also answer the central 
question about the origin of high-temperature superconductivity. 

The recent discovery of quantum oscillations in YBa2Cu306+a; (YBCOg+a;) with 
X ~ 0.5 and in YBa2Cu408 indicates that the Fermi liquid theory may also be applicable 
to underdoped cuprates O [71 El 13 HH], but the small size of the observed pockets 
is inconsistent with the large hole-like Fermi surface predicted by ab-initio electronic 
structure calculations. This suggests a reconstruction of the Fermi surface by an 
electronic superstructure, at least in the presence of the high magnetic fields required for 
the quantum oscillation experiments. An incommensurate magnetic superstructure has 
indeed been observed in underdoped La2_sSr^Cu04 (LSCO) [HI [T21 [131 [Il|, but until 
recently had appeared to be a peculiar property of this family of compounds. Moreover, 
quantum oscillations have not yet been observed in LSCO, presumably due to a higher 
level of intrinsic disorder. 

We recently reported a related superstructure in a YBCOe+x crystal with x = 0.45 
[T^ . slightly below the doping level at which the quantum oscillations were observed. 
Although the incommensurability is about a factor of two smaller than the one in LSCO 
at the same doping level, this implies that incommensurate magnetic order should be 
regarded as another generic zero-temperature phase in the cuprate phase diagram. It 
remains to be established whether the incommensurate Bragg reflections observed by 
neutron scattering are due to a modulation of the amplitude ("stripes") [16j or the 
direction ("spirals") [17] of the ordered moments, or both. Nonetheless, we refer to 
this state as a "spin density wave" (SDW), taking account of the observation that 
the incommensurate magnetic order in YBCO6.45 develops out of a metallic state, in 
contrast to the commensurate antiferromagnetic order in the Mott-insulating state at 
lower doping levels. The amplitude of the SDW modulation is strongly enhanced by a 
modest magnetic field H ~ 15 T [18]. It is thus very likely that the same SDW state 
is present at the much higher fields used in the quantum oscillation experiments on 
crystals with slightly higher doping levels. A reconstruction of the Fermi surface by the 
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SDW is then expected on general grounds. Subsequent work has confirmed that the 
size of the Fermi surface pockets observed in YBCOg.s can be quantitatively explained 
based on the incommensurate wave vector we have determined [I9l [20] . 

An investigation of the temperature dependence of the magnetic structure and 
dynamics of YBCO6.45 has revealed that the magnetic order disappears in a two-step 
fashion upon heating [15]. In a first step, muon spin relaxation (/xSR) experiments 
showed that static magnetic order, which reduces the translational symmetry of the 
lattice, vanishes at a temperature T ~ 2 K. In a second step, the incommensurate 
wave vector characterizing low-energy magnetic fluctuations continuously approaches 
the commensurate antiferromagnetic ordering wave vector with increasing temperature, 
and the incommensurability vanishes in an order-parameter- like fashion at Telc ~ 150 
K. While the magnetic fluctuations above Telc respect the fourfold rotation symmetry 
of the Cu02 plane, their wave vector spontaneously aligns with one of the two nearly 
equivalent copper-oxygen bond directions upon cooling below this temperature. In a 
wide temperature regime, therefore, the spin system collectively behaves in such a way 
that the rotational symmetry is broken while the translational symmetry remains intact. 
Since the symmetry properties of this state are analogous to those of a nematic liquid 
crystal, we refer to this state as an "electronic liquid-crystal" (ELC) [21]. The small 
orthorhombicity of the YBCO6.45 crystal lattice is required to align the ELC domains 
such that the ELC state becomes observable in a volume-integrating method such as 
neutron scattering. In the presence of this aligning field, the ELC transition is rounded 
into a sharp crossover. While the neutron scattering data are not accurate enough to 
quantify this rounding, the strong, order-parameter-like temperature dependence of the 
incommensurability indicates that the loss of rotational symmetry reflects an underlying 
phase transition that is driven by collective interactions between the spins. 

Our work on YBCO6.45 raises several questions. How is the incommensurate SDW 
related to the magnetic order observed in LSCO, and how does it evolve into the 
commensurate antiferromagnetic order [TJ [22] observed in YBCOe+z at lower doping 
levels? How is the ELC state related to anomalies in the charge dynamics observed 
by transport (such as Nernst effect measurements [231 [21]) and spectroscopic probes 
(such as infrared spectroscopy [2S]) in the same temperature and doping range? This 
article describes neutron scattering measurements on untwinned YBCOe+x crystals with 
X = 0.3 and 0.35 designed to address these questions, and to constrain theoretical models 
of the collective magnetic ordering phenomena in the cuprates. In conjunction with 
measurements on x = 0.45 and higher doping levels, the results allow us to determine 
the doping dependence of the incommensurate wave vector as well as the characteristic 
crossover temperatures for SDW and ELC ordering in YBCOe+x- For a YBCOe.ss 
crystal, we also present data on the magnetic field dependence of the SDW order 
parameter at low temperatures, as well as temperature dependent neutron resonant spin- 
echo experiments that bridge the gap in energy scales between /xSR and conventional 
neutron diffraction data. The results of our study allow us to draw an outline of the 
phase diagram of underdoped YBCOe+x as a function of doping, magnetic field, and 
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temperature. 

2. Experimental details 

The measurements were performed on crystal arrays with three different hole doping 
levels in the strongly underdoped regime of the YBCOg+x phase diagram. The samples 
are listed in Table [T} The oxygen content of the single crystals was carefully adjusted by a 
thermal treatment during which oxygen diffusion takes place. A variation of temperature 
at fixed oxygen partial pressure results in different oxygen contents and thus different 
hole concentrations [26]. The oxygen content and hole concentration per planar Cu 
ion, p, were extracted from the known doping dependence of the out-of-plane lattice 
parameter c and Tc [27j . 



Sample 


T 

^ c 


lattice parameter c ( A) 


hole doping p 


total mass mtot (g) 


YBC06.45 


35 K 


11.761 


0.082 


2.0 


YBC06.35 


lOK 


11.781 


0.062 


1.4 


YBCOe.s 


OK 


11.791 


0.052 


1.0 



Table 1. Summary of the underdoped YBCOe+a; samples. The SC transition 
temperature Tc was determined by magnetometry. The out-of-plane lattice parameter 
c at room temperature was determined by X-ray powder diffraction on selected samples 
from the same batch. From this value the hole doping level per planar Cu ion p could be 
extracted [17j- Several tens of crystals were co-aligned by X-ray Laue diffractometry, 
mimicking a large single-domain crystal with a total mass mtot- 

In YBCOg+x, a unique axis in the Cu02 plane is established by the orthorhombic 
distortion of the crystal structure. However, as-grown samples show crystallographic 
twinning, because the difference of the in-plane lattice parameters a and h is small. 
In order to observe the anisotropy of the electron system with a volume- averaging 
probe such as neutron scattering, a single-domain state has to be prepared in which 
the orientation of the two in-plane axes is maintained throughout the entire sample 
volume. This was accomplished by individually detwinning single crystals of YBCOg+a: 
(with typical size ax6xc = 2x2x 0.5 mm^) by application of uniaxial mechanical 
stress along the crystallographic (1,0,0) direction. 

The crystals were characterized by magnetometry and were found to exhibit SC 
transitions with widths of 2-4 K, testifying to their high quality. In order to obtain 
a sample volume sufficient for inelastic neutron scattering experiments, several tens of 
crystals were co-aligned on a single-crystalline silicon disc by X-ray Laue diffractometry. 
The resulting crystal arrays, mimicking large single-domain crystals, had mosaicities of 
< 1.5° and a majority twin domain population close to 90%. The contribution of the 
minority twin domains to the observed magnetic scattering signal is thus below the 
detection limit. 
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The neutron experiments were performed at the triple-axis spectrometers INS and 
IN14 (ILL, Grenoble, France), 4F1 and 4F2 (LLB, Saclay, France), Rita-II (SINQ, PSI 
Villigen, Switzerland), and on the resonant spin-echo triple-axis spectrometer TRISP 
(FRM-II, Garching, Germany). Pyrolytic-graphite crystals were used to monochromate 
and analyze the neutron beam. In order to extinguish higher-order contaminations of 
the neutron beam, a beryllium filter was inserted into the beam for those measurements 
which were performed with fixed final wave vector kf = 1.55 A . For the same 
purpose a pyrolytic-graphite filter was used when measuring at kf = 2.57 A and 
kf = 2.66 A , respectively. On TRISP, two radio-frequency coil sets, located between 
the monochromator and sample and between sample and analyzer, were used to 
manipulate the Larmor phase of the neutron spin before and after the scattering event, 
and the relaxation rate of magnetic excitations was extracted from spin-echo profiles 
[28]. 

The scattering vector Q = {H, K, L) is expressed in reciprocal lattice units (r.l.u.), 
i.e. in units of the reciprocal lattice vectors a*, b*, and c* {a* = 27r/a, etc.). As in 
the undoped (x = 0) parent compound, magnetic intensity is concentrated around the 
antiferromagnetic wave vector Qafm = (0.5, 0.5, L). Scans were thus performed either 
around Qafm around the equivalent points (0.5, L5,L) and (L5,0.5,L) in higher 
Brillouin zones. 

For the measurements in high magnetic field the sample was mounted in a 15-T 
vertical-field cryo magnet. The scattering plane was spanned by the vectors (1,0,0) and 
(0,1,2), and scans were performed along a* around Q = (0.5,0.5,1). The scattering 
geometry implies an angle of 33° between the magnetic field H and the c-axis. Thus, 
for the field-dependent experiments, the external magnetic field has a major component 
perpendicular to the Cu02 planes. All scans for H 7^ were performed after field 
cooling. 

3. Results 

3.1. Doping dependence 

We first discuss the behavior of the spin system at the lowest temperatures studied in 
our experiments (T = 2 K). Figure [l] shows the low-temperature neutron scattering 
function S{Q,u) of YBCO6.35 around the antiferromagnetic ordering wave vector 
Qafm = (0-5, 0.5, 1) as a function of excitation energy E = hco. It is apparent that 
the spectrum is comprised of two components: a sharp, intense "quasi-elastic" peak 
centered aX E = 0, and a weaker inelastic component that depends smoothly on energy. 
The two components are analogous to the Bragg peak and the spin wave spectrum of an 
antiferromagnetically ordered state. At the lowest temperatures, the energy width of the 
quasielastic peak is at the resolution limit of neutron resonant spin-echo spectroscopy 
(~ IfieV, see below), and /iSR experiments indicate magnetic order that persists much 
longer than the muon lifetime (~ 2/isec). At temperatures of order 2 K, which includes 
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Figure 1. Magnetic neutron scattering function of YBCOe.ss at T = 2 K and Q = 
(0.5, 0.5, 1), showing distinct quasielastic and inelastic components. 

the temperature range in which the quantum oscillation experiments were carried out, 
the quasielastic peak can therefore be regarded as a signature of static magnetic order. 
The two-component form of the magnetic spectrum is common to all three samples 
studied her^ 

In order to determine the spatial character of the magnetic correlations, we have 
performed constant-energy scans through the low-energy spin excitation spectrum. 
Figure [2] shows a comparison of constant-energy scans for E = ?> meV along the in- 
plane crystallographic directions a* and h* for YBCO6.45, YBCOe.ss, and YBCOg.s. 
Note that these are raw data, without any background subtraction. Scans through the 
quasielastic peak are of lower quality because they are superposed by a background 
due to incoherent scattering from the sample and the sample mount, but the scan 
profiles are nearly identical [I5]. Clearly, the magnetic response of all three samples 
shows a pronounced in-plane anisotropy. The cuts through the YBCO6.45 spectrum 
(Fig. |2^,b) demonstrate that the anisotropic intensity distribution is a consequence 
of two incommensurate peaks that are symmetrically displaced from Qafm along a*, 
whereas the distribution is commensurate along h* . For YBCOg.ss the anisotropy is 
qualitatively similar, but less pronounced (Fig. [2|:;,d), and it is further reduced for 
YBCOg.a (Fig. [2^,f). Note that even in this latter case we can rule out any influence 
of the instrumental resolution on the anisotropic peak profil^ Although the splitting 

I For YBCOe.a and YBCOg.ss, both commensurate quasielastic reflections centered at integer L and 
incommensurate peaks with continuous L-dependence characteristic of two-dimensional magnetic order 
are detected. Both sets of reflections exhibit a similar temperature dependence. 

§ The fuU-width-at-half-maximum of the elastic resolution function was determined at Qafm — 
(0.5,0.5,1) using the A/2 harmonic of the Bragg reflection (1 1 2) by removing the filters; it is 
~ (0.025 ± 0.002) r.l.u. in the a*- and - (0.023 ± 0.002) r.l.u. in the 6*-direction. The effective 
inelastic resolution is somewhat larger, but still much smaller than the width of the scan profiles both 
along a* and b* . 
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Figure 2. Constant-energy scans (raw data, without background correction) at i? = 
3meV and T = 3K along a* and b* for (a,b) YBCO6.45, (c,d) YBCOg.ss, and (e,f) 
YBCOe.s- The hues are single Gaussian (&*) and double-Gaussian (a*) profiles fitted to 
the data. The two incommensurate peaks symmetrically displaced from Qafm along 
a* are shown as dashed lines. 



between the two incommensurate peaks in the a* direction is no longer resolved for 
YBCOe.ss and YBCOg.s, fits to two-Gaussian profiles with individual peak widths 
matching the one in the b* direction provide excellent descriptions of the data (lines 
in Fig. [2]). Together with the observation of a split profile for YBCO6.45 (Fig. |2^), this 
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Figure 3. Doping dependence of the incommensurability 6 (i.e. the displacement 
of the peaks from Qafm) fo^' energy transfer of £^ = 3meV at low temperatures, 
extracted from the fits in Fig. |2] (closed symbols). The open symbol denotes the 
incommensurability of the low-energy spin excitations for Zn-substituted YBCOg.e 
(from [29|)- For comparison, the linear (5-versus-p relation for LSCO is shown as a 
dashed line \m\m. 



indicates an incommensurability of the magnetic response that increases continuously 
with increasing hole content in the Cu02 planes. 

Figure [3] summarizes the doping dependence of the incommensurability S in the 
three YBCOe+x samples that exhibit static magnetic order. We have also included 
data obtained from a YBCOg.e sample in which static incommensurate magnetic order 
was induced by substituting 2% of the Cu atoms by spinless Zn impurities [29]. It is 
instructive to compare these data to the doping dependence of 6 in the LSCO family 
of cuprates, where the propagation vector of the incommensurate modulation of the 
spin system is along the Cu-0 bond direction for p > 0.05, as it is in YBCOg+x [TTl [T2] . 
However, 6 is systematically lower in YBCOe+x than it is in LSCO at the same p. While 
in LSCO a linear relationship between the incommensurability and the doping level 
(5 = p) is observed for p < 0.12, S in YBCOe+x extrapolates to zero as p — )■ 0.04 from 
above, close to the doping level at which the propagation vector switches from the Cu-0 
bond direction to the diagonal of a CUO4 plaquette in LSCO [131 [E]. In contrast, we 
have not found any evidence for a 45° reorientation of the incommensurate modulation 
in YBCOg+x in the doping range we have investigated. Possible origins of this materials- 
specific behavior include exchange interactions between directly adjacent copper oxygen 
planes in the crystallographic unit cell (which is comprised of bilayer units in YBCOe+x 
and a single-layer unit in LSCO), and the lower disorder in YBCOg+a; which allows 
metallic conduction at lower p than in LSCO. Calculations considering the impact of 
these factors in the framework of a model in which the magnetic incommensurability 
arises from spiral order are indeed consistent with the behavior we have observed |30j . 
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Figure 4. Magnetic field dependence of the elastic magnetic peak at T = 2K for 
YBC06.45 (reproduced from ^8j) and YBCOe.as- The data were normalized to the 
intensity at H = 0. 



3.2. Magnetic field dependence 

As discussed in Section 1, we have shown that the intensity of the quasielastic peak 
in YBC06.45 is strongly enhanced in an external magnetic field, while its shape is only 
weakly affected [IB]. This effect is analogous to a similar phenomenology previously 
estabhshed in LSCO for a; ~ 1/8 [311 ESI |33j. It is naturally explained as a consequence 
of a competition between SDW and ci-wave superconducting phases, which coexist in 
YBC06.45 for H = [M]- Since the magnetic field destabilizes d-wave superconductivity 
by orbital depairing, the SC order parameter is reduced, and the SDW order parameter 
consequently enhanced for H 7^ 0. 

In order to further explore the validity of this scenario, we have extended our 
high-magnetic-field measurements to YBCOg.ss- Fig. |4] shows the H-dependence of the 
quasielastic peak intensity at T = 2 K for the two different doping levels. Whereas the 
peak intensity for YBCO6.45 increases linearly with increasing magnetic field [TH|, the 
magnetic intensity is field-independent within the error bar for YBCOe.ss- The field- 
independence is consistent with earlier observations on a twinned YBCOe+x sample with 
similar hole concentration [53] and on LSCO samples in which the SDW order is already 
well established for H = [36j. It also fits well into a two-phase competition scenario 
in which the SDW phase fraction in YBCOe.ss is already close to 100% for H = 0. 

3.3. Temperature dependence 

3.3.1. Spin density wave We now address the temperature dependence of the spin 
correlations, beginning with the intensity of the quasielastic peak. Figure |5^ shows the 
temperature dependence of the quasielastic peak intensity of YBCOe.ss measured on a 
conventional cold-neutron triple- axis spectrometer with energy resolution ~ 0.1 meV 



Neutron scattering study of the magnetic phase diagram of underdoped YBa2Cu'iO(i+x^^ 



YBCOg35 Q=(0.5, 0.5, 1), E=OmeV 




10 20 30 40 
Temperature T{K) 



5 10 15 20 25 30 35 
Temperature T{K.) 



Figure 5. Intensity of tlie quasielastic peak; in YBCOg.ss determined by conventional 
cold-neutron triple-axis spectrometry (a) and energy width (HWHM) of the same peak 
determined by resonant spin-echo spectroscopy (b) as a function of temperature. The 
lines are guides-to-the-eye. The inset shows a sketch of the Lorentzian lineshape of the 
quasielastic peak determined by spin-echo spectroscopy at different temperatures (solid 
lines) in relation to the Gaussian resolution of the triple-axis spectrometer (dashed 
line) . 



(half- width-at-half- maximum, HWHM). The data appear to show a SDW transition 
at a temperature of Tsdw ~ 40 K. However, the rounding of the transition, which 
substantially exceeds the width of the SC transition in these samples (~ 2-4 K), and 
the lack of saturation at low temperatures are unusual. Moreover, fiSK experiments 
indicate "spin freezing" transitions at much lower temperatures in this range of the 
phase diagram [13 [221 EZ]- Similar observations have been reported for underdoped 
LSCO [381 EH]. 

In order to explore the origin of this behavior, we have carried out neutron resonant 
spin-echo measurements [28j with an energy resolution of ~ 1 /xeV, which bridge the 
gap in energy scales between conventional neutron diffraction and /xSR. The spin- 
echo profiles indicate that the quasielastic peak exhibits a Lorentzian lineshape with 
an intrinsic energy width (HWHM), F, that exceeds the resolution of the spin-echo 
instrument for T > 5 K (Fig. |5}d). F increases continuously with increasing temperature 
and exceeds the HWHM of the Gaussian resolution function of the triple-axis instrument 
for T > 40 K (inset in Fig. |5^). These results show that the rounding observed in Fig. 
|5^ is of dynamic origin, and not due to an inhomogeneous distribution of transition 
temperatures as one might have suspected. They also imply that the SDW transition 
suggested by Fig. |5^ is actually a crossover whose specific form is determined by the 
energy window sampled by the experimental method. For comparison with the results 
of transport experiments and other experimental probes of the charge dynamics, it is 
nonetheless useful to monitor the strength of the quasi-static SDW correlations at T 7^ 0, 
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which can be characterized by the crossover temperature Tsdw above which the width of 
the quasielastic peak exceeds the energy scale ~ 0.1 meV. We have therefore extracted 
this temperature from the data of Fig. |5^ and analogous data collected for the other 
doping levels. 

The diverging fluctuation time scale implied by the data of Fig. [5]d is naturally 
explained as a consequence of the zero-temperature phase transition expected for a 
two-dimensional spin system with full spin-rotation symmetry. In this scenario, the 
intensity of the quasielastic peak for T — )• is proportional to the SDW order parameter. 
In contrast to the divergence of the correlation length expected for a T = phase 
transition, the width of the scattering profiles in momentum space (Fig. |2]) saturates at 
low temperature. This may reflect either an inhomogeneous broadening of the profiles 
due to slight variation of doping levels across the large single-crystal array, or finite-size 
effects imposed by the random distribution of dopant atoms and/or the coexistence 
of SDW and SC phases. We note, however, that the two-component form of the 
scattering function shown in Fig. [T] is closely similar to the spectrum of an undoped 
two-dimensional square-lattice antiferromagnet at T 7^ 0, where thermal fluctuations 
generate a finite, intrinsic correlation length [101 E]. Possible origins for the intensity 
decrease of the inelastic component as — )■ in the low-temperature data of Fig. [T] 
include disorder-induced finite-size correlations and the intra-atomic spin-orbit coupling, 
which leads to anisotropy gaps in the spin- wave spectra of undoped cuprates p]. 

3.3.2. Electronic liquid crystal At temperatures exceeding Tsdw ~ 40 K, the 
quasielastic peak of Fig. [l] is no longer visible, but prominent low-energy 
incommensurate spin fluctuations persist. Since the incommensurability could not be 
resolved in most of the inelastic scattering profiles for x < 0.35, we have extracted 
the momentum widths (FWHM) of the profiles in the a* and h* directions, and 
Afo, from scans like those shown in Fig. [2] at different temperatures. Whenever the 
incommensurate peak separation 25 exceeded A;,, A^ was extracted from a two-Gaussian 
fit. Otherwise the a* scan was fitted by a single Gaussian for the sake of simplicity. For 
either way of fitting, A^ is the difference of the H values for which the intensity drops 
to half of the maximum intensity. 

Figure |6] shows the resulting widths Aq and A;,, corrected for the instrumental 
resolution. While A^ increases gradually upon heating, A^ decreases and continuously 
approaches A;,. The temperature dependence of Aq was therefore fitted by a polynomial 
function, while for A^ a simple linear fit was applied which accounts for the thermal 
broadening of the peaks. The point of intersection of these two curves marks Telc? 
the onset temperature of the incommensurability. Above Telc? both widths are nearly 
identical and increase in a parallel manner upon further heating. 

In order to compare the three different samples directly, we have plotted the 
temperature dependence of the anisotropy, defined as the deviation from unity of 
the ratio of A^ and the linear fit of A^ (Fig. [?]). While the low-temperature limit 
of the anisotropy decreases with decreasing hole concentration, reflecting the doping 
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Figure 6. Temperature evolution of the effective momentum widths (FWHM) along 
the a and b directions of the spin excitation profiles at i? = 3meV for (a) YBCO6.45 
(b) YBCO6.35, and (c) YBCOe.a- The widths were extracted from constant-energy 
scans, as shown in Fig. [2] and corrected for the instrumental resolution. 
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Figure 7. Comparison of the temperature-dependent anisotropy, Aa/Af, — 1, of the 
effective momentum widths A^ and A;, for E — 3mcV for the three different hole 
doping levels. The temperature for which the polynomial fit yields Aa/A;, — 1 = is 
defined as the ELC transition temperature Telc- 



dependence of the incommensurability (Fig. [2]), the temperature where Aa/A;, — 1 — t- 
increases with decreasing p. 

The strong temperature dependence of the in-plane anisotropy of the spin 
fluctuations in all three samples and its appearance at a well-defined, doping dependent 
onset temperature confirms that this phenomenon arises from collective interactions 
between spins. The "Ising nematic" transition expected on general grounds for a two- 
dimensional incommensurate magnet [211 1121 SSI HH SS], possibly associated with a 
"Pomeranchuk" instability of the Fermi surface ^61 W7\ HE], provides an interesting 
qualitative explanation of our observations, but a quantitative assessment of this 
scenario will have to await further theoretical progress. Meanwhile, it is interesting 
to compare our YBCOe-i-a; data to analogous observations on the Ba2Fe2-a;(Co,Ni)a;As2 
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system, where a strong, doping dependent in-plane anisotropy of the spin excitations 
has also been found |l9l |50]. In contrast to the cuprates, however, this anisotropy 
is temperature independent and can be associated with nesting features of the Fermi 
surface [50]. The SDW and ELC states in YBCOg+a;, which are observed in immediate 
vicinity of the Mott insulator, thus do not appear to be straightforwardly related to the 
more conventional SDW state in the iron pnictides. 

4. Summary 

Figure [8^ summarizes the overall layout of the p-T phase diagram of underdoped 
YBCOg+x- The SDW and ELC crossover lines established by our neutron scattering 
results are shown along with the previously determined phase boundaries for 
commensurate antiferromagnetism [22] and superconductivity [SJJ- At T = 0, the 
SDW and SC phases coexist over some range of p. Our data are consistent with /zSR 
evidence of coexistence between magnetic order and superconductivity in this doping 
regime [221 ES] , but /xSR does not contain information about the spatial character of the 
local magnetization. The neutron scattering data presented here demonstrate that the 
magnetic order coexisting with superconductivity is incommensurate. At the present 
time, they are insufficient to determine whether the coexistence is spatially uniform, as 
predicted by spiral models [JTIEO], or spatially modulated in the form of stripes [TUIIHH] 
or spin vortex lattices [M]- In any case, unavoidable oxygen defects in the YBCOe+x 
crystal lattice are expected to induce some degree of spatial inhomogeneity. 

Both the SDW order parameter extracted from the intensity of the quasielastic 
peak as T — i- and the crossover temperature Tsdw that characterizes its onset at 
T 7^ are continuously reduced with increasing p. This indicates that the SDW phase 
boundary ends at a quantum critical point (QCP) at p ~ 0.1, close to the chemical 
composition YBCOg.s. Further evidence for quantum criticality is derived from the 
scaling properties of the dynamical spin correlations of YBCO6.45 [55]. In line with this 
scenario, the scattering function at T = exhibits a spin gap for p > 0.1 [56], but a 
dilute concentration of nonmagnetic impurities restores the SDW |29] . 

The low-p end of the incommensurate SDW regime is complicated by a confluence 
of several factors including the tetragonal-to-orthorhombic transition [2B] , the onset of 
appreciable effective interactions between CUO2 bilayers (~ where J±_ ~ 0.02 

meV is the bare inter-bilayer interaction [1] and ^ is the intra-layer correlation length, 
which increases with decreasing p), and dopant-induced disorder, which may lead to 
the formation of local Cu^"*" moments on the CuO chains for samples with low oxygen 
content [FT] . Nonetheless, it is interesting to note that the onset of incommensurability 
in the spin correlations (Fig. [3]) nearly coincides with the onset of superconductivity 
(Fig. [s]). The YBCOe+x system is less dominated by disorder than strongly underdoped 
LSCO ^58] and thus provides an interesting complementary forum for future exploration 
of the interplay between spin and charge correlations close to the Mott insulating state. 

In order to relate our low-temperature neutron scattering data to the results of the 
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Figure 8. Magnetic pliase diagram of underdoped YBCOg+a;. Panel (a) shows a p-T 
section of the phase diagram with the crossover temperatures Tsdw and Telc extracted 
from Figs. [5^ (and similar data for other samples) and[7| respectively, along with the 
phase boundaries of commensurate antiferromagnetism [221 and superconductivity |51) . 
Panel (b) shows a schematic of a p-H section of the phase diagram derived in a field- 
theoretical treatment of competing SDW and SC order parameters [34]. The control 
parameter in the field theory was heuristically identified with the doping level p. The 
field range covered by the data of Fig. |4] and the larger field range covered by the 
quantum oscillation experiments on YBCOe.s samples are indicated by red arrows. 



quantum oscillation experiments, we have schematically drawn the range of magnetic 
fields covered by our experiments on top of the zero-temperature phase diagram for 
coupled SDW and d-wave SC parameters derived from field theory ^] (Fig. [8|d). 
Here, the control parameter in the field theory was heuristically identified with the hole 
concentration p, without claims to theoretical rigor or quantitative accuracy. Whereas 
YBCOg.ss appears to be near the SDW end of the coexistence regime, far away from 
the QCP, both our YBCO6.45 sample and the YBCOe.s samples used in the quantum 
oscillation experiments are close to the QCP. The limited magnetic fields currently 
available for neutron scattering are incapable of promoting the electron system into 
the uniform SDW phase expected at high H, but sufficient to substantially increase 
the SDW order parameter in the coexistence regime. In models with uniform SDW- 
SC coexistence for H = 0, this proceeds via a spatially periodic enhancement of the 
SDW order in the vicinity of magnetic vortices in the mixed state of the d-wave SC. 
It will be interesting to explore a possible relation between the spatially modulated 
SDW order and the series of field-induced transitions recently observed in the magnetic 
vortex lattice of optimally doped YBCOe+x [59]. The much higher fields used for the 
quantum oscillation experiments are expected to bring the electron system of samples 
with doping levels near the zero-field QCP either close to the high-field boundary of the 
SDW-SC coexistence range or into the uniform SDW phase (Fig. [sjo). An extrapolation 
of our low-field measurements of the SDW order parameter of YBCO6.45 to H = 50 T 
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yields ~ O.lSfiB per Cu ion [18], comparable to the zero-field sublattice magnetization 
of stripe-ordered Nd- or Eu-substituted LSCO, where recent photoemission experiments 
have uncovered evidence of a reconstructed Fermi surface [60l l6T] . It is thus likely that 
a SDW reconstruction of the multiband Fermi surface of YBCOe.s contributes to the 
formation of the small pockets observed in the quantum oscillation experiments. This 
conclusion is supported by recent quantitative assessments of the quantum oscillation 
data [iniEn]- It is interesting to note that the QCP indicated by our neutron scattering 
study coincides with the insulator-metal transition inferred from transport experiments 
in high magnetic fields [621 [63], and that the transport data indicate a critical divergence 
of the cyclotron mass at this doping level [531 El] • A comprehensive explanation of the 
neutron-scattering and quantum-oscillation data therefore requires a transition from 
"SDW insulator" to "SDW metal" at high fields [US], for which there is thus far no 
direct experimental evidence. 

Going back to the zero-field phase diagram of Fig. [8^, we note that the ELC 
onset temperature also increases with decreasing p, such that the ELC and SDW 
crossover lines are nearly parallel. In view of the minute size (~ 0.5% for YBCOg.ss) and 
opposite doping trend of the orthorhombic distortion of the crystal structure [2SIES], this 
observation supports our contention [15] that the anisotropy is not a simple consequence 
of the crystal symmetry [HZ], but arises spontaneously from many-body interactions 
between spins. The ELC onset temperature is substantially lower than the onset of 
the "charge pseudogap" detected by infrared spectroscopy [251 EHl [69], which at these 
doping levels is well above room temperature. However, it coincides with the onset of 
signatures of phase-incoherent superconductivity in the infrared spectra of YBCOe+i 
samples from the same batch as ours ^D]. The theoretical description of the subtle 
interplay between superconducting and spin correlations for T ~ Telc is an interesting 
subject of further investigation, as is the relationship between the ELC phenomenon 
reported here and the spontaneous onset of in-plane anisotropics in do resistivity [66] 
and Nernst effect [23] measurements. While the temperature evolution of the in-plane 
resistivity anisotropy [SS] is in good qualitative agreement with the anisotropy of the 
spin dynamics shown in Fig. ^ corresponding Nernst-effect data are not yet available 
for the doping levels investigated here. 

Finally, it will be important to follow the ELC crossover line to higher doping levels, 
where the SDW phase is no longer present at T = 0. For YBCOe.e {p ~ 0.12), the spin 
excitation spectrum exhibits a gap in the superconducting state at T = [56j. Weak 
low-energy spin excitations are present above Tc, but their temperature dependence has 
not yet been studied in detail. The intensity of higher-energy excitations with E 30 
meV, however, shows a strong upturn upon cooling below ~ 150 K [7^, in qualitative 
agreement with transport [HSl 123] and infrared [HH] data. Taken together, these results 
suggest that the ELC phenomenon persists at least up to p ~ 0.12. It will be interesting 
to explore the relationship between this transition and the "spin pseudogap" detected 
by nuclear magnetic resonance [7T], as well as the magnetic refiections with Q = 
detected by elastic neutron scattering [721 EH] in this doping regime. 
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